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Abstract 

This study is an investigation of spin wave spectrum in macrostructures composed of two fer- 
romagnetic materials and showing a 3D periodicity: spherical ferromagnetic grains disposed in 
the nodes of a 3D crystal lattice are embedded in a matrix with different ferromagnetic properties. 
Frequency ranges forbidden to spin wave propagation are found in the calculated magnonic spectra. 
Both the position and the width of the gaps are found to depend on the magnetic (exchange and 
magnetization) contrasts in the composite material, as well as on its structural parameters (filling 
fraction, crystal lattice type). Having applied our theory to interpretation of the existence of a spin 
wave gap in doped manganites, recently revealed in neutron scattering experiments by S. Hennion 
et al., we obtained a good (though approximate) quantitative agreement with the experimental 
results. A working hypothesis is proposed on this basis, supposing that certain manganites can be 
treated as natural 3D magnonic crystals. 
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INTRODUCTION 



Though the first study of electromagnetic wave propagation in periodic structures, writ- 
ten by lord Rayleigh, was published already in 1887, it is only in recent years that photonic 
composites raised suddenly an extremely keen interest. The research in this field was ini- 
tiated by the studies by Yablonovitch and John published in 1987 and anticipating 
the existence of complete energy gaps in electromagnetic wave spectra in three-dimensional 
periodic composites, henceforth referred to as photonic crystals. These can be used for 
fabricating new optoelectronic devices in which the role of electrons, traditionally used as 
transport medium, would be played by photons Q,[ [ 0] . The so-called left-handed materials 
(LHM), showing negative effective refractivity |a, llf, provide an example of periodic mate- 
rials demonstrating how much the properties of this kind of structure can differ from those 
of homogeneous materials. Another type of periodic composites are structures composed 
of materials with different elastic properties; showing an energy gap in their elastic wave 
spectrum, such composites are referred to as phononic crystals Recently, attention 
has been focused on the search of photonic and phononic crystals in which both the position 
and the width of the energy gap could be controlled by external factors, such as applied 
voltage or magnetic field. Attempts are made to create photonic crystals in which one of 
the component materials would be a magnetic 

A magnetic periodic composite consists of at least two magnetic materials; the informa- 
tion carrier in such structures are spin waves. By analogy to photonic and phononic crystals, 
in which the role of information carrier is played by photons and phonons respectively, pe- 
riodic magnetic composites are referred to as magnonic crystals. Studies of 2D magnonic 
crystals have been reported with scattering centers in the form of "infinitely" long 

cylinders disposed in square lattice nodes (cylinder and matrix materials being two different 
ferromagnetics), and the anticipated gaps found indeed in the respective spin wave spectra. 
Further research was focused on magnetic multilayer systems, which can be regarded as ID 
magnonic crystals 

In this paper, we present numerically calculated magnonic band structures - to our best 
knowledge not yet reported in the literature - of three-dimensional magnonic crystals. Due 
to the complexity of the problem, only the simplest model of 3D magnonic crystal is consid- 
ered here, represented by a system of ferromagnetic spheres (which act as scattering centers) 
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disposed in the nodes of a cubic crystal lattice embedded in a different ferromagnetic ma- 
terial (matrix). Both the exchange and dipolar interactions are taken into account in our 
calculations, based on the plane wave method and using the linear approximation. As a con- 
clusion, we propose a new magnonic interpretation of ex per imental results obtained through 
neutron scattering on spin waves in doped manganites j^^j-j^l- 



THE PROCEDURE OF NUMERICAL CALCULATION OF 3D MAGNONIC BAND 
STRUCTURE 



Let's consider an ideal periodic structure consisting of spheres of ferromagnetic A and 
embedded in a matrix of ferromagnetic B. The spheres are assumed to form a 3D periodical 
lattice of sc (Fig. [T^) or bcc (Fig. [TJo) type. A static magnetic field, Hq, is applied to the 
composite along the X3 axis, and assumed to be strong enough to saturate the magnetization 
of both materials. The lattice constant is denoted by a; the filling fraction, /, is defined as 
the volume proportion of material A in a unit cell: 

f = ^ (sc) or f = ^ (bcc), 1 
in an sc or bcc lattice, respectively. 



a) sc b) bcc 




FIG. 1: The 3D periodic structures studied in this paper; two types of structure are considered, both 
consisting of ferromagnetic spheres of material A embedded in a matrix of material B (materials 
A and B have different magnetic properties); (a) the spheres are disposed in the nodes of an sc 
lattice; (b) the spheres are disposed in the nodes of a bcc lattice. 



Ferromagnetics A and B are characterized by two material parameters: the spontaneous 
magnetization {Ms^a and Ms^b), and the exchange constant {A^ and Ab)] both these pa- 
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rameters depend on the position vector r = (xi,X2,X3): 

M5(r) = Ms,B + {Ms,A - Ms,B)S{f), 
A{r) = AB + {AA-AB)S{r), (2) 

the value of function S{r) being 1 inside a sphere, and beyond. 

In the classical approximation, spin waves are described by the Landau-Lifshitz (LL) 
equation, taking the following form in the case of magnetic composites (with damping ne- 
glected) : 

dMir t) ^ ^ 

^ = 7/UoM(f,t) X H,ff{r,t), (3) 

where magnetization M{f, t) is a function of position vector r and time t; 
Heff{f, t) stands for the effective magnetic field acting on magnetization M(r, t) and defined 
as follows 



i?e//(r, t) = Hoe, + K(f, + ^ (v ■ M(f, t); (4) 

Cj is the unit vector along the Xj axis; /I(r, t) is the dynamic magnetic field due to dipolar 
interactions; the third component represents the exchange field. The magnetization vector 
can be represented as the sum of its static component, Mse,-, which is parallel to the applied 
field, and its dynamic component, m{f,t), lying in plane {xi,X2)'- 

M{f, t) = Mse3 + m{f, t). (5) 

The dynamic dipolar field, /i, must satisfy the magnetostatic Maxwell equations: 

V X /I(r) = 0, 

V ■ (h{r) + m(f)) = 0. (6) 



In magnonic crystals, the position-dependent coefficients in Q, i.e. Ms and A, are periodic 
functions of the position vector, which allows us to use, in the procedure of solving the 
Aj equation defined in 0, the plane wave method, described in detail in our earlier papers 
(dealing with two- dmensiona/ magnonic crystals). Following this scheme, we proceed 
to Fourier-expanding all the periodic functions of the position vector, i.e. the spontaneous 
magnetization, M^, and parameter Q defined as follows: 



The dynamic components of the magnetization can be expressed as the product of the 
periodic envelope function and the Bloch factor exp{iKf) {K denoting a 3D wave vector); 
the envelope function can be transformed into the reciprocal space as well. Including all 
the expansions into and (jH]) leads to the following infinite system of linear equations for 
Fourier coefficients of the dynamic magnetization components, rh^^{g) and m2i:{g)'- 

-no ^ \k + q'r 0. 

9' 9" ^ ' 

l-K It: 
■ Ms{-{g - g"))Q{-{g" - ^))m,0), 

yjj (Jj 

•o , ^ \^ (M + g^}(M + gi) ,27t 

-no ^ \k + g'\^ a 
-no ^ \k + g'\^ 0- 



[9 -9") ■(9-9')] 



9' 9" 

■ Ms{-{g-g"))Q{-{9"-9'))m,d9')- (8) 

(Jj (Jj 

ki and gi denoting the Cartesian components of a dimensionless wave vector k = {a/27c)K 
and a dimensionless reciprocal lattice vector g = {a/2n)G, respectively; a new quantity 
introduced in (jH)) is Q, henceforth referred to as reduced frequency: 

n = (9) 

The Fourier coefficients of spontaneous magnetization Ms and parameter Q are calculated 
from the inverse Fourier transformation; in the case of spheres the resulting formulae read 
as follows: 

' Ms,Af + Ms,B{^-f), for G = 

Ms{G) = \ , (10) 
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and similarly 



Q{G) 



QAf + QB{l-f), 



for ^ = 



(11) 



where Qa and Qb are the values of parameter Q in material A and B, respectively; R is 
the sphere radius (Fig. [TJ. 

Obviously, the numerical calculations performed on the basis of (jH)) involve a finite number 
of reciprocal lattice vectors G in the Fourier expansions; however, we have made sure that 
the number used is large enough to ensure good convergence of the numerical results. As 
indicated by an analysis performed, a satisfactory convergence is obtained already with 343 
reciprocal lattice vectors used. 

a) sc b) bcc 





FIG. 2: . The Brillouin zones corresponding to the considered structures, and the paths (high- 
Ughted) along which the magnonic band structures are calculated: (a) MTXMR in the case of the 
sc lattice, (b) PTHNP in the case of the bcc lattice. 



NUMERICAL RESULTS 

The first 3D magnonic crystal to be studied here is a magnetic composite consisting of 
ferromagnetic spheres (of material A) disposed in the nodes of a simple cubic [sc] lattice 
and embedded in a different magnetic material (B), referred to as matrix (Fig. ^). The 
corresponding magnonic band structure will be calculated along path M = 7r/a(l, 1,0) — >■ 
F = 7r/a(0,0,0) -> X = 7r/a(l,0,0) M = 7r/a(l,l,0) R = 7r/a(l,l,l) in the 

6 



nonreducible part of the first Brillouin zone (see Fig. ^.). Iron and yttrium iron garnet 
(YIG) are chosen as component materials A (spheres) and B (matrix), respectively, in the 
studied example. As established in our earlier studies jl?!. in the case of 2D magnonic 
crystals, such composition, involving a substantial contrast of the magnetic parameters 
between the component materials [Ms^Pe = 1.752 lO^Am"^, Mgyic = 0.194 lO^Am"^, 
Apf. = 2.1 10~^^ Jm~^ , Ayig = 0.4 10^^^ Jm~^), leads to large energy gaps in the spin wave 
spectrum. Also in the spin wave spectrum obtained for the sc lattice-based 3D magnonic 
crystal considered here, two wide energy gaps are present, one between bands 1 and 2, the 
other between bands 5 and 6, as shown in Fig. IHl The existence of these gaps signifies that 
spin waves having frequency values within a gap cannot propagate in the composite. The 
following parameter values have been assumed in our calculations: crystal lattice constant 
a=100A, applied static magnetic field /io-f^o = O.IT, and filling fraction /=0.2. 



sc lattice 

400 1 

G 




I ' 

M r X M R 

Wave vector, k 



FIG. 3: The magnonic band structure found numerically for the sc lattice-based composite (iron 
spheres embedded in a YIG matrix). The spin wave energy branches have been calculated along the 
path (in the first Brillouin zone) shown in Fig. 12^. The following parameter values are assumed: 
applied static magnetic field ^jlqHq = O.IT, crystal lattice constant a = lOOA, and filling fraction 
/=0.2. 

Let's examine the effect of the filling fraction on the width of the gaps obtained. Figure 
I3K shows the magnonic bands plotted against the filling fraction in the sc lattice-based 
structure (Fe spheres embedded in YIG); the plot was obtained through projecting the 
magnonic band structure calculated along MTXMR (for a fixed filling fraction value) on 
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the reduced frequency (Q) axis; the procedure of projection was performed repeatedly for 
consecutive filhng fraction values. The first (lowest) gap shown in Fig. is found to be the 
largest for filling fraction value /=0.24, the corresponding gap width being Ar2=63.29, while 
the much narrower second gap reaches its maximum width at /=0.3. Unlike in ID and 2D 
magnonic crystals, no oscillatory variations of the energy gaps with the filling fraction are 
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22j | . The gaps found in the bcc lattice-based structure, 



observed in the 3D composites 
shown in Fig. |3Jd, are clearly much broader than those in the sc lattice-based structure. The 
spin wave reduced frequency was calculated along path PTHNP [P = 7r/a(l,l,l) — > F = 
7r/a(0,0,0) H = 7r/a(2,0,0) ^ = 7r/a(l,l,0) ^ P = 7r/a(l,l,l)] (indicated by a 
bold line in Fig. |2b) in the bcc lattice first Brillouin zone. The maximum width of the 
first gap is found to occur at /=0.2, reaching value Ai7=l 17.02, which is 85% larger than 
its counterpart in the sc lattice-based structure. The second energy gap found in the bcc 
lattice-based structure (between bands 4 and 5) is broader than the first one, and reaches a 
maximum width of AQ =141.52 at filling fraction value /=0.32. 

a) sc lattice b) bcc lattice 
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FIG. 4: The magnonic energy branches in cubic magnonic crystals (Fe spheres embedded in YIG) 
plotted versus the filling fraction (a) for an sc lattice-based structure, and (b) for a bcc lattice-based 
structure. The lattice constant value is assumed to be a = lOOA, and the applied magnetic field 
value to be HoHq = O.IT. 



Let's proceed now to the role played by the magnetic parameters of the component 
materials. Figure El shows spin wave energy branches plotted against the contrast between 
the spontaneous magnetization values in materials A (spheres) and B (matrix); this contrast, 
defined as ratio Ms^a/Ms^b, will be henceforth shortly referred to as magnetization contrast. 
The computations have been performed for a fictitious matrix material whose spontaneous 



8 



magnetization and exchange constant values, M^ ^ and Ab, respectively, are fixed at values 
close to those in YIG; the exchange constant value in the sphere material is assumed to be 
Aa = 2.1 lO^^^Jm^^. The other parameters are fixed as well: the assumed lattice constant 
value is a=100A, the filling fraction /=0.2, and the applied magnetic field fioHo = O.IT. The 
results obtained for the sc lattice-based and the bcc lattice-based structures are depicted in 
Figs. 5a and 5b, respectively. Let's focus on the first gap occurring between the first and the 
second energy bands; the gap width is found to be maximal for magnetization contrast values 
10 and 12.8 in the bcc lattice-based and the sc lattice-based structures, respectively. Both 
these values are not far from the spontaneous magnetization contrast between iron and YIG 
{Ms^Fe/ MsyiG = 10.82), the material composition used in the above-discussed investigation 
of the effect of the filling fraction, and leading to quite wide gaps, as one can see in Figs. 
and EJd. Note that in both lattice types the minimum spontaneous magnetization contrast 
necessary for energy gaps to open is greater than 2. 



a) sc lattice b) bcc lattice 




0.0 4.0 8.0 12.0 16.0 0.0 4.0 8.0 12.0 16.0 

Magnetization contrast, Magnetization contrast, 

Ms,a/ Ms.b Ms.a/ Ms.b 



FIG. 5: Three-dimensional magnonic crystal spin wave spectra plotted versus the magnetization 
contrast, on the basis of computations performed for (a) an sc lattice-based structure, and (b) a bcc 
lattice-based structure. The assumed spontaneous magnetization and exchange constant values in 
the matrix are Ms^b = 0.194 W^Am~^ and Ab = 0.4 10~^^ Jm"^, close to the respective parameter 
values in YIG; the exchange constant value in the spheres is assumed to be Aa = 2.1 10^^^ Jm^^, 
the lattice constant a = lOOA, the filling fraction / = 0.2, and the applied field ^qHq = O.IT. 

In the last stage of this investigation, we shall examine the role of the contrast between the 
exchange constant values in the component materials. The results of the respective computa- 
tions are depicted in Fig. showing magnonic branches versus the exchange contrast, defined 
as the ratio [Aa/Ab) between the exchange constant values in materials A (spheres) and B 
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(matrix); results obtained for sc and hcc lattice-based structures are shown in Figs. IHK and 
inb, respectively. The only variable in the computations was the exchange constant in mate- 
rial A, all the other structural and material parameters being fixed: Ms^a = 1-752 10^ Am'^, 
Ms,B = 0.194 W^Am-\ Ab = 0.4 10^^^Jm'\ a = lOOA, ^o^^o = O.IT. The occurrence 
of the first energy gap (between the first and the second magnonic bands) is found to be 
independent of the exchange contrast, yet the gap width grows as the exchange contrast 
increases. This means that the exchange contrast is not an indispensable factor for the 
opening of complete energy gaps in 3D magnonic crystals, its effect being limited to the 
width of the existing gaps. 

a) sc lattice b) bcc lattice 
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FIG. 6: Magnonic branches plotted versus the exchange contrast in the composite. The assumed 
values of the matrix material parameters are close to the respective parameter values in YIG; the 
spontaneous magnetization in the spheres is assumed to be Ms^a = 2.1 10^ Am~^. Other parameter 
values are as in Fig. 5. 



CONFRONTATION OF THE THEORY WITH THE EXPERIMENT 

As demonstrated by the above-presented results, complete energy gaps can occur in the 
spin wave spectra of 3D periodic magnetic composites. The opening of such gaps is condi- 
tioned on the material composition, as gaps are found to exist only when the spontaneous 
magnetization contrast is greater than 2. The gap width can be controlled through adjusting 
the filling fraction value, and the position of the gap on the energy scale through modifying 
the value of the applied magnetic field. As expected, the bcc lattice-based structure is found 
to give rise to larger energy gaps than the sc lattice-based one. 
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FIG. 7: Magnonic band structure computed on the basis of Eqs © for a hypothetic periodic 
structure of ferromagnetic droplets existing in doped manganites Lai-xCuxMnO^ (1% < x < 
9%); the droplets are assumed to be sphere-shaped and to form a regular crystal lattice of bcc 
type. The material and structural parameter values used in the computations are as follows: 
Ms,A = 0.8 106Am-^ Ms,b = 0.2 lO^Am-^ Aa = 0.04 10^lVm-^ Ab = 0.02 l^-^^Jm-^, 
f = 0.25, a = 50 A, ^qHq = O.IT, close to the respective estimate values obtained experimentally 



We shall try to confront our theory with available experimental results now. Energy 
gaps have recently been found in spin wave spectra of doped manganites Lai -T.Cacr.MnO?. 
and Lai-xSr^MnO^ at low concentrations (2% < x < 9%) of Ca and Sr ions S^-js^. 
The quoted studies demonstrated the existence of ferromagnetic nanoregions (referred to as 
droplets) in these materials at temperatures below a critical point; the space between the 
droplets is filled with a different magnetic (ferromagnetic or antiferromagnetic) medium. We 
suggest that the spin wave spectrum gap recently found experimentally through neutron scat- 
tering on doped manganites might be an indirect indication of periodic spatial distribution 
of the ferromagnetic droplets. Assuming that the droplets (material A) are sphere-shaped 
and disposed in the nodes of a bcc lattice, we have computed the magnonic spectrum in 
a 3D composite with structural and material parameter values close to those reported in 
the experimental studies The resulting spin wave spectrum is shown in Fig. d 

the following parameter values are assumed: Ms^a = 0.8 10^ Am~^, Ms^b = 0.2 10^ Am~^, 
Aa = 0.04 10'" Jm-i, Ab = 0.02 10-"Jm-\ / = 0.25, a = 50A, /igi^o = O.IT. A narrow 
energy gap is found in the obtained spectrum (between the first and the second bands), in 
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ihe same frequency range as the gap which was found in the neutron scattering experiments 
S^-js^. The gap resulting from our numerical computations lies between 1.27meV and 
1.49meV; the gap found experimentally had width 0.2meV and was reported to lie in the 
vicinity of 1.65meV. As the "computed" gap has width 0.22meV, the result obtained on 
the ground of our theory can be deemed close to the experimental one. This preliminary 
result fills us with optimism, as it allows us to propose a working hypothesis supposing that 
the manganites in question can be regarded as 3D manganite structures. This, however, 
will remain unconfirmed as long as no direct experimental evidence of the regularity of the 
structure formed by the droplets is available (c/. results reported in j^sl)- We shall continue 
our research in order to find further evidence supporting our hypothesis. 
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